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Neutrino oscillograms of the Earth and CP violation in neutrino 
oscillations 
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69117 Heidelberg, Germany 

An analysis of 3-fiavour neutrinos oscillations inside the Earth is presented in terms of the oscillograms - 
contour plots of oscillation probabilities in the plane neutrino energy - nadir angle. Special attention is paid to 
CP violation in neutrino oscillations in the Earth. 



1. INTRODUCTION 

Neutrino oscillograms of the Earth are contours 
of constant oscillation or survival probabilities in 
the plane of neutrino energy and nadir or zenith 
angle of neutrino trajectory. They encode infor- 
mation on both the neutrino parameters and the 
Earth density profile and proved to be a very use- 
ful and illuminating tool for analyzing neutrino 
oscillations in the Earth. 

The oscillograms exhibit a very rich structure 
with local and global maxima and minima (in- 
cluding the MSW resonance maxima in the man- 
tle and core of the Earth), saddle points and the 
parametric enhancement ridges (Fig. []}. It was 
shown in p] that all these features can be read- 
ily understood in terms of various realizations of 
just two conditions: the generalized amplitude 
(collinearity) and phase conditions. The study 




in [T] was performed in the limit Am^ = 0; 
in the present talk, based on [2], the results for 
Am|x 7^ are presented. 



2. OSCILLOGRAMS FOR 3-FLAVOUR 
NEUTRINO OSCILLATIONS 

In the full 3-flavour framework with Amfj ^ 
and 0i3 7^ 0, the oscillograms have non-trivial 
structure both at low and high energies. Os- 
cillations at low energies are dominated by the 
solar channel, whereas those at higher energies 
are mainly driven by the atmospheric parameters 
(An&.tfis) (seeFigrj. 



Figure 1. Oscillograms for 1 — P ee for three different 
values of #13. 



The main qualitatively new feature of 3-flavour 
oscillations as compared to 2-flavour ones is the 
possibility of CP violation. Its effect can be con- 
veniently studied in terms of CP-oscillograms, 
which are defined as contour plots of equal prob- 
ability difference 



AP, 



CP 



P»e(6) - P^ e (S th ) . 



(1) 



where S is the true value of the CP violating phase 
and 8th is the assumed, or theoretical, value which 
we want to test. The CP-oscillograms have a 
rather complex domain- like structure (see Fig. [J); 
however, this structure can be readily understood 
and interpreted in terms of three grids of curves, 
as we discuss next. 
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Figure 2. Oscillograms for the difference of prob- 
abilities APC P = P Me (5)-iV(<5th) with J th = 0°. 
Shown are the solar (black), atmospheric (white) 
and interference phase condition (cyan) curves. 



The oscillation probability P Me can be conve- 
niently written as 

P^e = \c 23 A e2 e* 5 + s 23 A e3 \\ (2) 

where the amplitudes A e2 and A e3 are indepen- 
dent of #23 and CP-violating phase S. To leading 
order in small parameters Si 3 and Ato^i/ Ato^, 
A e2 depends only on the solar oscillation param- 
eters and A e3 , only on the atmospheric ones: 

A e2 ~ A s (6 12 , Am 2 21 ) , A e3 ~ A A (6 13 , Am§ x ) , 

(3) 

i.e., these amplitudes approximately coincide 
with, respectively, the solar and atmospheric con- 



tributions to the amplitude of «-> oscilla- 
tions. From eq. © we find 

= c 2 23 \A s \ 2 + s 2 23 \A A \ 2 

+2s 23 c 23 \A s \\A A \cos(<p + 5), (4) 

where = arg(AgA A ). Only the last term in (T?)) 
depends on 5; the condition of vanishing AP^ e 
(i.e. P Aie ( ( 5) = P^e(<5th)) therefore corresponds to 

\A s \\A a \cob(4>+S) = \A s \\A A \coB(4>+St h ). (5) 

There are three non-trivial realisations of this 
condition: 

• As = (=> solar "magic" lines) 

• A A = atm. "magic" curves) 

• (<t> + Sth) = -{<t> + S) + 2Trl, or 
ct> = -(5 + 5 th )/2 + nk. 

There is also, of course, the trivial realization, 
Sth = S + 27m, when the true and assumed values 
of S coincide. The above three conditions give 
rise to three grids of curves, which delineate the 
different domains in the CP oscillograms (Fig. [2]) ; 
they have a very simple physical interpretation. 
We shall now discuss these conditions in the en- 
ergy region between the 1-2 and 1-3 MSW reso- 
nances; the results for other energy domains can 
be found in [2]. 

2.1. Solar "magic" curves 

The condition of vanishing solar channel con- 
tribution to the oscillation probability, As = 0, 
in the constant matter density approximation can 
be written as sin 2 (<x>i2£) = (i.e. uj± 2 L = 7m), 
where lui 2 is the solar oscillation frequency. At 
energies exceeding the 1-2 resonance energies in 
the mantle and in the core of the Earth, E > 0.5 
GeV, one has lo 12 ~ V/2, where V = y/2G F N e 
is the matter-induced potential of neutrinos. The 
condition As = therefore takes the form 



where L(Q V ) is the nadir angle dependent length 
of the neutrino trajectory in the Earth. Note 
that the condition ^ is energy independent and 
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determines the baselines for which the "solar" 
contribution to the probability vanishes. In the 
plane (Q V ,E V ) it represents nearly vertical lines 
G„ w const (black lines in Fig. [2]). 

There are three solar magic lines which cor- 
respond to n = 1 (in the mantle domain) and 
n = 2,3 (in the core domain). The existence of 
a baseline (L w 7600 km) for which the proba- 
bility of v e <-» oscillations in the Earth is ap- 
proximately independent of the "solar" parame- 
ters (Amjji, #12) and of the CP-phase 8 was first 
pointed out in [3] and later discussed in a number 
of publications (see |4|5j and references therein) . 
This baseline was dubbed "magic" in [4|. The 
interpretation of this baseline as corresponding 
to vanishing "solar" amplitude A e 2 , according to 
eq. ((6|) with n — 1, was given in [5]. It was also 
shown there that for neutrino trajectories cross- 
ing the core of the Earth there exist two more 
solar "magic" baselines, corresponding to the os- 
cillation phase equal irn with n = 2 and 3, and 
the existence of the atmospheric "magic" curves 
was pointed out. 

2.2. Atmospheric "magic curves" 

The second realization of the condition |5]) 
corresponds to vanishing atmospheric amplitude, 
Aa = 0. In this case, like on the solar "magic" 
lines, the probabilities of v e <-> (as well as 
z/ e <-> j/ r ) oscillations are independent of the CP- 
violating phase. In addition, they do not depend 
on the atmospheric parameters Am^ and $13. 

The properties of atmospheric "magic" curves 
can be easily understood in the constant den- 
sity approximation, in which the condition Aa = 
is satisfied when sin 2 (a;3iL(6;,)) = 0, i.e. 

W3iL(8„) = irk, k = 1, 2, For energies which 

are not too close to the atmospheric MSW reso- 
nance energy, this condition reduces to 



E 



Am§ 1 £(9 i/ ) 
\Airk±2VL(Q v )\ 



(7) 



which corresponds to the bent curves in the 
(€>!,, E v ) plane (shown in white in Fig. [2]). 

2.3. The interference phase condition 

The third of the above mentioned realizations 
of the condition ^ depends on the interference 



phase between the amplitudes Ag and Aa and 
therefore we shall call it the interference phase 
condition. In the energy region between the 1-2 
and 1-3 resonances we have 4> ~ Am^L/AE, i.e. 
in the first approximation <fr does not depend on 
the matter density. The interference phase con- 
dition then takes the form Am^L/AE = — (6 + 
6 th )/2 + irl, or 



E v = 



Am| 1 L(e y ) 



(8) 



Anl - 2{S + 6th) ' 

From the comparison of eqs. ([5]) and (J7J) it follows 
that the interference phase curves (shown in cyan 
in the oscillogram of Fig. \2§ are similar to atmo- 
spheric curves, but are steeper than the latter for 
v's and less steep for D's. 

3. CONCLUSIONS 

It can be seen from Fig. [2] that the three grids 
of curves discussed above - the solar and atmo- 
spheric "magic" ones and the interference phase 
curves - give a very accurate description of the 
borders between the domains. On these bor- 
ders the probability difference AP^ e p = -P Me (<5) — 
-P^ie(^th) vanishes, and therefore there is no sen- 
sitivity to the phase <5. The regions of maximal 
sensitivity to the CP-violating phase correspond 
to the central parts of the domains. 

Our analysis has been performed in terms of 
the oscillation probabilities; a more realistic study 
which considers the number of events in future 
detectors is currently under way. 
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